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Abstract: The detection of the flavour content of jets produced from R-parity 
violating neutralino decays is investigated in the case where one baryon-number 
violating coupling dominates. Simulations are performed of the ATLAS experiment 
at the LHC for all couplings, other than A" jfc since neutralino decays through these 
couplings are very suppressed. Secondary vertex distributions and muons produced 
by heavy-quark (b- and c-) jets allow discrimination between LSP decay modes. The 
dominant coupling can be identified at better than 3.5 o in almost all cases, with the 
only remaining ambiguity caused by the inability to distinguish strange from down 
quarks. 
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1. Introduction 

In supersymmetry (SUSY) models where R-parity, R p = (— i) 3B + L + 2S ; i s violated 
(RPV models) the lightest supersymmetric particle (LSP) can decay to Standard 
Model particles. For the baryon-number violating RPV coupling, A", the LSP, as- 
sumed to be the Xi, decays via a virtual squark into three quarks (or antiquarks). 
In previous work [1] we suggested a method by which the masses of several of the 
SUSY particles can be measured in such scenarios for a wide range of A". 

If R-parity violating couplings are observed, it seems probable that they come 
from an underlying theory which also determines the fermion mass spectrum. This 
idea has been explored [2-5] in the context of gauged U(l)x models of flavour, in 
which a hierarchical pattern of RPV couplings is generated by the same symmetry 
from which the fermion mass spectrum is derived. Thus the determination of the 
flavour structure of an RPV coupling could be a vital clue in the resolution of the 
flavour problem. 

In this paper we investigate the flavour content of the quark jets from Xi decays. 
The aim is to identify which X"j k coupling is dominant in the neutralino decay. The 
index i is the generation number of an up-type quark, and j and k correspond to 
down-type quarks. The coupling is antisymmetric in j and k, so there are nine 
possible non-zero elements. 

Three of the couplings, A^ fe , produce very different signatures because they allow 
the Xi t° decay via the top quark. If the neutralino is lighter (heavier) than the top 
quark then the decay is propagator (phase-space) suppressed. For the mSUGRA 
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Figure 1: The lifetime of the xi at the mSUGRA point described in the text, with one 
of \'-j k non-zero. The solid line shows the lifetime for a dominant top quark coupling, 
which leads to the decay Xi ~ > qq't* — > qq'bW* — ► qq'bff' , where q and q' are quarks other 
than top and f and f are Standard Model fermions. The dashed line shows the lifetime 
for couplings to quarks other than the top. The calculation approximates the final state 
particles as massless, so that the lines for \"- k are the same for all j, k, as are those for 
for all i / t, j, k. 

point we investigate (which is described in Section 0) the branching ratio through 
the top coupling is a factor of 4.2 x 10 6 smaller than through a non-top coupling of 
the same size. 

For a top mode to dominate the decay, a A" coupling to top would need to be 
a factor of at least 2000 greater than any of the non-top couplings. In that case 
the lightest neutralino would typically be long-lived (Figure [l|) and for A" jfe ^ 10~ 2 
nearly all Xi s decay outside the detector. Determination of these couplings could 
then only be achieved by searching for the rare decay of Xi s m the active volume, 
while sparticle mass measurements could be made in the same way as for the R p 
conserving case [6]. 

We investigate the other six couplings - A" ds , A" db , A" sb , A" ds , A" db and A" sb - 
assuming that one is significantly larger than the others. Since the branching ratio 
is proportional to the square of the coupling, it is sufficient for this analysis that all 
the subdominant couplings are 4.5 times smaller than the dominant coupling. Then 
just 5% of all x? s will decay though the subdominant channel, and 90% of events, 
each with two x?s, will contain only dominant-coupling decays. 
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2. Event Simulation and Selection 

HERWIG 6.2 [7-10] is used as the event generator, and the official ATLAS simulation 
program ATLFAST 2.50 [11], is used to simulate the performance of the ATLAS 
detector. We use as a test case the mSUGRA point m = 100 GeV, m\j% = 300 GeV, 
A Q = 300 GeV, tan/? = 10 and sgn/i +, with one of Af jfc = 5 x 10~ 3 . This is in the 
middle of the range of couplings 10~ 5 < A" < 1, for which we can expect double x? 
decay within the detector, and is approaching indirect limits set from nuclear decays 
for Ay ds and A^ db [12-14]. A selection of the sparticle masses for this point is shown in 
Table [l|. The applicability of our analysis to SUGRA models with other parameters 
is discussed in our previous work [1]. The method of the event selection, which is 
also detailed in that paper, is summarised below. 

Signal events will contain three jets from each of the two Xi s - Since the SUSY 
cross section is dominated by squark and gluino production, there are usually at least 
two other jets from the R p conserving decays of squarks and gluinos. Requiring the 
presence of leptons, which can originate in the decay chain x® ~> ^r^ ± ~~ * Xi-^^ > 
helps to decrease the Standard Model (SM) background. Events are selected with 
between eight and ten jets, and with two opposite-sign, same family isolated leptons 
(e or /i). 

Cuts on event-shape variables are applied, reducing the SM background to less 
than 10% of the SUSY signal, and kinematic cuts are applied to preferentially select 
jets from neutralino decays. All possible three-jet combinations are inspected and 
their invariant masses calculated. Since we expect the event to contain two LSP 
decays we select two of these three-jet combinations with reconstructed masses within 
20 GeV of each other as Xi candidates. To reconstruct the x\ we Pi°k the Xi candidate 
closest in r]—(j) to either of the leptons from the lepton pair, and find the total invariant 
mass of the three jets plus those two leptons. 

A 2-dimensional Gaussian is fitted to the Xi an d x\ masses, and combinations are 
selected within 2 x a of that peak. This selection removes much of the combinatoric 
background, so that the jet multiplicity cut can be relaxed to 8 < iVj Ct < 12. We 
now have a rather clean sample of events in which we can identify six jets with the 
two neutralino decays. All iV com b jet combinations which pass the above cuts for any 
event are accepted and are given weight l/iV comb . 



X? 


X°2 


g 


fii? 




d R 


~d L 


Ir 


h 


116.7 


211.9 


706.3 


611.7 


632.6 


610.6 


637.5 


155.3 


230.5 



Table 1: Masses of selected particles (in GeV) for the model investigated. These were 
calculated from the mSUGRA parameters using ISASUGRA from the package IS A JET 7.51 
[15]. 
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Figure 2: The fraction of signal events with 0, 1, and 2 vertex tags, for different types of 
RPV coupling A". The statistical errors are those expected for an integrated luminosity of 
30 fb _1 . Points are horizontally displaced to allow them to be distinguished by eye. 



3. Vertex Tagging 

The lifetimes of hadrons containing b- and c-quarks give rise to displaced vertices 
which can be reconstructed from charged tracks in the inner detector. In general 
the lifetimes of charmed hadrons are shorter than for hadrons containing a b-quark, 
for example cr^o = 124 /zm compared to cr#o = 468 /im. This allows statistical 
separation of c-quark jets from b-quark jets, and of c- and b-jets from light quark 
(u, d, s) and gluon jets. Strange hadrons have longer lifetimes with cr of the order 
of tens of centimeters. Since they decay to a small number of particles it is difficult 
to reconstruct a secondary vertex, and so strange hadrons cannot be tagged within 
jets in the LHC environment. 

The vertex tagging performance of the ATLAS Inner Detector was simulated 
in [6]. In that study a likelihood ratio method was applied to the transverse impact 
parameter (do) of the tracks within the jet cone. The resultant b-tagging efficiency 
and the rejection rates for c- and other jets were parameterised as functions of the 
transverse momentum and pseudorapidity of the reconstructed jet and were imple- 
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b-tagging efficiency 


Rejection factor 


u, d, s and g jets (r,) 


c-jets (r c ) 


0.33 


1400 


22.9 



Table 2: b-tagging efficiencies and mis-tagging rejection factors from a full simulation 
of the ATLAS inner detector [6]. l/r c is the probability of tagging a c-jet, while 1/rj is 
the probability of tagging a u, d, s or gluon jet, averaged over pr, r\ and 4>. The typical 
reconstructed jet pr scale is 50 GeV. 

mented in ATLFAST. The rejection factors, averaged over all directions and trans- 
verse momenta, are summarised in Table In order to select a higher purity sample, 
we tolerate a rather low b-tagging efficiency (33%). 

The proportions of events with 0, 1 and 2 tagged jets are shown in Figure |2], for 
the six different A" couplings. As we expect, the vertex tagging rate is greatest when 
the X? decay products include both c and b quarks, and smallest for couplings which 
produce light-quark daughters. As we cannot distinguish d- from s-jets, there is no 
discrimination between A" db and X'- sh (for i=u,c). 

In these simulations, with A" = 5 x 10~ 3 the Xi lifetime of 10 -14 seconds corre- 
sponds to a typical decay length cr of 3 fim. This is short in comparison with typical 
b- and c-jet vertex displacements so it should not severely affect vertex reconstruc- 
tion. However, the lifetime is inversely proportional to the square of the coupling 
strength. When cr becomes bigger than about 3 cm (for A" ^ 5 x 10 -5 ) special 
vertex reconstruction would be required. Vertex tagging becomes almost impossible 
when neutralinos travel more than about 30 cm in the transverse direction since their 
daughters will produce hits only in the outermost silicon layers. 

4. Flavour Discrimination from Muons 

The weak decay of hadrons containing heavy quarks can produce muons, which 
will generally lie within the associated quark jet. The muons will pass though the 
calorimeter and so can be measured by the muon detector even if they lie inside 
the jet cone. The frequency with which these "non-isolated" muons occur in signal 
events can be used as an additional discriminator between quark flavours. 

The number of non-isolated muons per event in the six signal jets is plotted in 
Figure |3|. Since bottom mesons produce more muons [B(B° — > fi + z/ M X) = 10%] 
than charmed mesons [B(D° — > [i + X) = 7%] [16], we can use the number of 
muons in signal jets to statistically separate b- from c-jets. 

5. Statistical Significance 

The confidence with which we can identify the dominant coupling was explored for 
all pairs of couplings for both the vertex tagging and muon rates (Figures |2| and |3] 
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Figure 3: The fraction of signal events with 0, 1 and 2 muons that lie within 
y A?? 2 — A0 2 < 0.4 of the signal jets, for different types of RPV coupling, A". The 
statistical errors are those expected for an integrated luminosity of 30 fb _1 . Points are 
horizontally displaced to allow them to be distinguished by eye. 



respectively). The variable, 



X 



E 



i X j 



was calculated, where X{ and yi are the fractions of events with i muons (or vertex 
tags) for the test pair of couplings and the a 2 s are their variances assuming Poisson- 
distributed numbers of events. Bins where one or other distribution contained fewer 
than five events were excluded. 

The calculated \ 2 values for both the muon and vertex-tagging plots are shown 
in Table |||. If we use only the muon information then all couplings are distinguishable 
at 90% except for the ambiguity between d- and s-jets. In general the vertexing rate 
gives better discrimination, but would be difficult for a long-lived Xi f° r t ne reasons 
discussed in Section |3]. Combining the results from both analyses gives separation 
at better than 3.5 a in all cases, again with the exception that it is not possible to 
distinguish down from strange quarks. 
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Table 3: Chi-squared function, and number of degrees of freedom (d.f.) for the difference 
in distributions for pairs of RPV couplings (A") for an integrated luminosity of 30 fb" 1 . The 
contributions from vertexing and muon counting are shown separately. The probability (P) 
in the tail of the chi-squared distribution is given when P > 1%. The number of degrees 
of freedom is one less than the number of histogram bins in which both couplings contain 
at least five events. 

If the branching ratios of the neutralino through two or more couplings were 
comparable, then as the information in the distributions is rather degenerate, defini- 
tive identification of the couplings becomes difficult. However the method could then 
be used to constrain the possible values of the couplings. 

6. Conclusions 

If baryon-number violating couplings to quarks other than top are > 5 x 10~ 5 
then x?s will typically decay within the tracking volume of an LHC detector. For 
these couplings, our simulations show that displaced secondary vertices and muons 
from heavy quark decays allow statistical separation of b- from c-quark jets, and 
b- and c-jets from light-quark jets. A dominant RPV coupling can then be identified 
at better than 3.5 a, except for an ambiguity caused by the inability to distinguish 
strange from down quarks. 

We have demonstrated that if the MSSM breaks baryon number, not only can 
the LHC detect and measure the masses of sparticles [1], but it can even distinguish 
the flavour structure of the RPV coupling. 
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